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■ Abstract 

. We propose an effective description of squark interactions with charginos/neutra- 

linos. We recompute the strong corrections to squark partial decay widths, and 
compare the full one-loop computation with the effective description. The effective 
description includes the effective Yukawa couplings, and another logarithmic term 
which encodes the supersymmetry-breaking. The proposed effective couplings repro- 
duce correctly the radiative-corrected partial decay widths of the squark decays into 
charginos and neutralinos in all relevant regions of the parameter space. 



1 Introduction 

The Standard Model (SM) of the strong and electroweak interactions is the present para- 
digm of particle physics. Its validity has been tested to a level better than one per mille 
at particle accelerators [1]. Nevertheless, there are arguments against the SM being the 
fundamental model of particle interactions [2] , giving rise to the investigation of compet- 
ing alternative or extended models, which can be tested at high-energy colliders, such as 
the Large Hadron Collider (LHC) [3,4], or a 500 — 1000 GeV, e^e~ International Linear 
Collider (ILC) [5,6]. One of the most promising possibilities for physics beyond the SM 
is the incorporation of Supersymmetry (SUSY), which leads to a renormalizable field the- 
ory with precisely calculable predictions to be tested in present and future experiments. 
The simplest supersymmetric extension of the SM is the Minimal Supersymmetric Stan- 
dard Model (MSSM) [7-10]. Among the most important phenomcnological consequences 
of SUSY models, is the prediction of new particles. There is much excitement for the 
possibihty of discovering these new particles at the recently built LHC [11, 12], and their 
properties will need to be precisely measured to confirm (or refute) that they belong to 
a SUSY model. This last effort might be better suited for the ILC [5,6,13-15], currently 
being projected. This job needs the performance of precision measurements, but also of 
precision computations which are well suited for experimental comparisons. In the present 
work we will focus on the properties of the SUSY partners of the SM quarks - the squarks. 

Once produced, squarks will decay in a way dependent on the model parameters (see 
e.g. [16]). If gluinos (the fermionic SUSY partners of gluons) are light enough, squarks 
will mainly decay into gluinos and quarks (g qg) [17, 18], which proceeds trough a 
coupling constant of strong strength. If the mass difference among different squarks is 
large enough, some squarks can decay via a bosonic channel into an electroweak gauge 
boson and another squark (§„ — > 9^(2', W^)), and if Higgs bosons are light enough, also 
the scalar decay channels are available (g^ Qbi^^^ H"^)) [19-22], which can be 

dominant for third generation squarks due to the large Yukawa couplings. Otherwise, 
the main decay channels of squarks are their partial decays into charginos/neutralinos 
(the fermionic SUSY partners of the electroweak gauge and Higgs bosons) and quarks 
{q q'x)- Some of those channels are expected to be always open, given the large mass 
difference between quarks and squarks, and that the charginos/neutralinos are expected to 
be lighter than most of squarks in the majority of SUSY-brcaking models. In the few cases 
in which these channels are closed, the squarks will decay through flavour changing neutral 
channels [23-25], or through three- or four-body decay channels involving a non- resonant 
SUSY particle [26-31]. 

Here wc will concentrate on the squark decay channels involving charginos and neutra- 
linos. Their partial decay widths were computed some time ago, including the radiative 
corrections due to the strong (QCD) [32-34], and the electroweak (EW) [35-38] sectors 
of the theory. These radiative corrections are large in certain regions of the parameter 
space [37] , and their complicated expressions are not suitable for their introduction in the 
montc-carlo programs used for experimental analyses. In this work we present approxi- 
mations for the partial decay widths of squarks into charginos and neutralinos, including 
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the QCD corrections, and compare these approximations against the fixed-order one-loop 
corrected partial widths. 

In section |2] we introduce our notation and conventions for particles and couplings, and 
set up the numerical values that we will use in our analysis, section [3] presents the QCD 
one-loop computation of the partial decay widths and shows some numerical examples, 
in section H] we perform a renormalization group analysis of the partial decay widths, in 
section [5] we perform a numerical comparison of the one-loop and renormalization group 
computations, and finally section [6] shows our conclusions. 



2 Notation, conventions and numerical setup 

To describe the computation of the partial decay widths, we will follow the conventions of 
Ref. [39]. Throughout this work we will use a third-generation notation to describe quarks 
and squarks, but the analytic results and conclusions are completely general, and can be 
used for quarks- squarks of any generation. We will show numerical results only for third 
generation quarks/squarks (top t/stop t/bottom 6/sbottom b), since their decay widths 
are the ones that present the most interesting properties. 

We will study the partial decay widths of sfermions into fermions and charginos/neutra- 
linos, 

r(/-/'x). (1) 

We denote the two sfermion-mass eigenvalues by mj^ (a = 1,2), with m^^ < mj^. The 
sfermion-mixing angle ^/ is defined by the transformation relating the weak- interaction 
if a = II, /j?) and the mass eigenstate (/„ = /i, /2) sfermion bases: 

f _ d(/) f> . d(/) _ / cos% -sin% \ 

By this basis transformation, the sfermion mass matrix, 

_( Ml+m} + c,,{T, - Q si) Ml ruj Mf^ 

■^^/"i nifMj^ Ml+mj + Qcpsl^Ml f ' ^"^^ 



becomes diagonal: R'--^^ M.'^^R^-^^'^ = diag < m| , m| >. M? is the soft-SUSY-breaking mass 
parameter of the SU (2) l doublelo, whereas M? is the soft-SUSY-breaking mass parameter 

^ J R, 

of the singlet. T3 and Q are the usual third component of the isospin and the electric charge 
respectively, mj is the corresponding fermion mass, and sw is the sinus of the weak mixing 
ang lei The mixing parameters in the non-diagonal entries read 



M^^ = Ab- fitanp , M/^^ = - /i/tan/3 . 



^With — due to SU{2)l gauge invariance. 

^We abbreviate trigonometric functions by their initials, like sw = sin 9w, C2/3 = cos(2/3), tw = sw/cw, 
etc. 
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Ah^t are the trilinear soft-SUSY-breaking couplings, /i is the higgsino mass parameter, 
and tan/? is the ratio between the vacuum expectation values of the two Higgs doublets 
tan/3 = V2/V1. The input parameters in the sfermion sector are then: 

(M;^, M^^,Mj^, A, A, /i, tan /5) , (4) 

for each sfermion doublet. From them, we can derive the masses and mixing angles: 

i^b^^^h^^b) , {mi^,mi^,et) . (5) 

For the trilinear couplings, we require the approximate (necessary) condition 

Al<3 (ml + ml + Mfj + i?) , (6) 

where rriq is of the order of the average squark masses for q = t,b, to avoid colour-breaking 
minima in the MSSM Higgs potential [40-43]. 

Although the tree-level chargino (x^)-neutralino sector is well known, we give here 
a short description, in order to set our conventions. We start by constructing the following 
set of Weyl spinors: 

r+ ^ Hw+,H^) , 

r- = {-tW-,H-) , (7) 

The mass Lagrangian in this basis reads 



C 



M 



4r.r-) 





M 



r+ 



(ri,r2,r3,r4)A^° 



r2 

Vr4/ 



+ h.c. 



where we have defined 
M = 

= 



f M V2MwS/3 
\V2MwCp /i 

/ M' 

M 

Mzcpsw -Mzcpcw 

\-MzS13Sw Mzspcw 



MzCpSw -MzSpSw\ 
-MzCpCw MzSpCw 







/ 



(9) 



with M and M' the SU{2)l and U{1)y soft-SUSY-breaking gaugino masses. The four- 
component mass-eigenstate fields are related to the ones in ([7]) by 



Xi 



+T 



OT 
a 1 
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where U, V and N are in general complex matrices that diagonalize the mass-matrices (Q: 

U*MV^ = Md = diag(Mi,M2) (0 < Ml < M2) , , . 

N*M°N^ = M% = diag (MO, MgO, ikfO) (0 < M° < M° < < MO) . ^ ' 

Using this notation, the tree-level interaction Lagrangian between fermion-sfermion- 
(chargino or neutralino) reads [37] 

a=l,2 r 

^..Lf = -9f:Xr(Af^rPL + A^IlPR)f'- (H) 



"Xrfaf 

Here we have adopted a compact notation, where /' is either / or its SU{2) l partner for 
Xr being a neutralino or a chargino, respectively. Roman characters a,b . . . are reserved for 
sfermion indices and i,j, . . . for chargino indices; Greek indices a, (3, . . . denote neutralinos; 
Roman indices r,s . . . indicate either a chargino or a neutralino. For example, the top- 
squark interactions with charginos are obtained by replacing /—> t, Xr ^ X 

r = 1,2. The coupling matrices that encode the dynamics are given by 

A(t) _ p(*)t/* \ t?(*)T/* 
A-ii = -XbRalUi2, 

= -^{Ri!liK, + YLtwN:,) + V2x,R^^iN:, 

^-ai - -XtRal Vi2 , 

41 = -^{RS{N:2-yLtwK,)-V2x,R^iN:,) , 

A^-L = -^(-v^A,4?iV,3 + F^twi?SiV,i) , (12) 

with Yl and Y^^ the weak hypercharges of the left-handed SU (2) ^ doublet and right- 
handed singlet fermion, and At = rrit /{\/2 Mw sin P) and A;, = mb/{V^Mw cos P) are the 
Yukawa couplings normalized to the SU{2)l gauge coupling constant g. Note the following, 
each coupling is formed by two parts: the gaugino part, formed exclusively by gauge 
couplings, and the higgsino part, which contains factors of the quark masses, each of these 
parts will receive different kinds of corrections (see below). 

Using these definitions, the tree-level partial decay widths read 

IbiTm- ■I'' ■' 

J a 

K - - m),) [\AfS + \A^ll?) -^rrif, M.Re {a^H A^PS)] ,(13) 
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with A(x^, ?/^, z"^) = ^/[x'^ — {y — zy][x'^ — {y + zy]- 



2.1 Numerical setup 

For the numerical analysis and plots we will use fixed values for the SUSY parameters, and 
make plots by changing one parameter at a time. For the central values of the parameters 
we take: 

tan/3 = 5 , /i = 300 GeV , M = 200 GeV , Mj^ = 800 GeV , m-g = 3000 GeV , 

MsusY = Mj^ = 1000 GeV , At = = 2Mj^ + ^/ tan/5 = 1660 GeV , ^ ^ 

where we have introduced a parameter Msusy as a shortcut for all the SUSY mass pa- 
rameters which are not explicitly given. We use the GUT relation M' = b/St^^M for the 
bino mass parameter. For the SM parameters we use = 171.2 GeV, = 4.7 GeV, 
a^iMz) = 0.1172, sly = 0.221, Mz = 91.1875 GeV, 1/a = 137.035989. The renormaliza- 
tion scale Q is taken to be the physical mass of the decaying squark. The value of the 
trilinear couplings Ah^t is given by the algebraic expression, the given numerical value cor- 
responds to the default values of the other parameters, this numerical value will change in 
the plots, the chosen expression allows to show plots with a significant parameter variation 
avoiding colour-breaking- vacuum conditions ([6]). With these input parameters, the central 
values for the physical SUSY particle masses are: 

= (170.40, 337.50) GeV , 

M^o = (89.52, 172.28, 305.46, 338.58) GeV , 

= (802.05, 1000.30) GeV , 

rui = (720.55, 1084.25) GeV . (15) 

It is illustrative to identify the largest EW-basis component in each physical state. Of 
course, we have performed our computation using the full numerical mixing among the 
EW-basis and the physical-basis components, but this identification will help us to analyze 
the numerical results. The lightest squarks (ti,6i) are predominantly left-handed, the 
lightest chargino and neutralinos {XijXijX2) are predominantly gaugino-like, whereas the 
heaviest ones (x^jXsjX^) are predominantly higgsino-like. Of course, the parameters in 
eq. (fn|) are just and example for illustrative purposes, we have checked that our conclusions 
hold for a wide range of the parameter space. 



3 QCD Corrections 

Following this setup, we have computed the full one-loop QCD corrections to the squark 
partial decay widths into charginos and neutralinos ([1]). The renormalization prescriptions 
follow that of Ref. [37]. The QCD corrections include contributions from gluon loops, 
gluino loops, and gluon bremsstrahlung. The full one-loop corrections have been performed 
using the FeynArts/FormCalc/LoopTools packages [39,44-48]. We have used dimensional 
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reduction for the regularization of ultraviolet (UV) divergences, and a small gluon mass to 
regularize the infrared (IR) divergences. The three-body phase-space integration of the real 
gluon emission is performed analytically over the full energy range, and the dependence 
on the gluon mass is seen to cancel between the virtual and the real corrections. We have 
found full agreement with previous works [33,34,37,38], and will not repeat the full lengthy 
formulae here. The corrections are seen to be numerically large, specially in certain regions 
of the parameter space [37], specially those involving processes with a bottom-squark in 
the initial state, and in a regime of large tan (3 values. 

We follow the hints from Higgs-boson physics [49-52], and define effective Yukawa 
couplings which should encode the leading part of the corrections [50]: 

^eff. _ ^1^^' _ mb{Q) 
\ = 



Vl vi{l + Amb) 



.eff. _ mf^- _ mtiQ) 

where mq{Q) is the running quark mass and Aruq is the finite threshold correction. The 
SUSY-QCD contributions to Am^ are: 

^^SQCD ^ '^rngfitanP I{mc,mc,mg) , 
ovr 

^^SQCD ^ J(^^-^^^_^^_) ^ (17) 

where the function J(a, b, c) is the scalar three-point function at zero momentum transfer, 
and reads: 

ln(a2 /P) + Pc^ ln{byc^) + c^a^ ln{c^ jo?) 

The effective description of squark decays consists in replacing the tree-level quark masses 
in the couplings ffT^ by the effective Yukawa couplings of eq. ffTBl) . and use this lagrangian 
to compute the partial decay width, schematically: 

YYuk.-eff. ^ rtree^^e//.^) _ (^g) 

This expression contains the large one-loop corrections from the finite threshold correc- 
tions (fT6|) . but it also contains higher order corrections. At this point we can make the 
following: make a computation that combines the higher order effects (which ignore the 
effects of external momenta) and the fixed one-loop (which ignore the higher order effects). 
At the same time, this will allow us to quantify the degree of accuracy obtained by the 
effective description. We define a Yukawa-improved decay width computation: 
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where 



-pi— loop 
^1—loop 

cYuk.—rem. 



-ptree _j_ ^-pl—loop 



ptree 

pi— loop pi— loop Yuk.—eff. 



'^tree i 



m 



eff. 



(20) 



Here r^~'°°^ is the one-loop fixed order prediction for the partial decay width, Y^~'-°°'^ Yuk.-eff. 
is the one-loop expansion of the prediction using effective couplings, and therefore, the re- 
mainder contribution ((J^^^ -^e™- ) jg the part of the one- loop contribution that can not be 
described by the Yukawa effective couplings, it quantifies the approximation done by the 
effective description. 

The one-loop effective prediction pi-'oop jg computed by taking the computa- 

tion using effective couplings f|T8l) . expanding it in series, and keeping only the one- loop 
terms. Specifically: 



m{Q) 



m(m) 



l--«.(Q)log(^) + ... 
TT m 



and 



m 



eff. 



m{Q){l - Am) 



therefore, the part of the one-loop effective mass is: 



m(m] 



2 Q 
— tts(Q) log(— ) - Am 
TT m 



(21) 



f22) 



(23) 



this is the mass that will be used in the effective Yukawa couplings to compute r^"''"'^' Yuk.-eff. 
and m{m) is the running quark mass at the quark mass scale. Finally, we define a Yukawa- 
improved correction factor in the following way: 



yYuk.—imp. 'ptree ( ^ \ 

^Yuk.-imp. _ ^ V"'q) 



r 



tree i 



(24) 



All these definitions will allow us to precisely analyze the approximations. As an 
example. Fig. [T] shows the partial decay width (and the relative correction) of a top- 
squark decaying into the lightest chargino, as a function of the SU{2)l squark mass scale 



([3]), the rest of the parameters are given in flMj) . We see a big dip in the corrections 
for squark masses around 1250 GeV, with negative corrections surpassing —100% - which 
would mean a negative decay width, which obviously does not make sense. What happens 
is that, for this very special setup of parameters, the tree-level computation of the partial 
decay width vanishes, so the one-loop contribution exceeds the tree-level prediction. Under 
these circumstances one-loop perturbation theory does not hold, and we can not claim the 
validity of any result obtained by the one-loop perturbative expansion, that is: we can not 
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(a) (b) 

Figure 1: a) Partial decay width and b) relative corrections to the top-squark decay width 
into the lightest chargino as a function of the SU{2)l squark mass scale Mj^, for the 
parameters of eq. f|T^ . 



give a prediction for the decay width in those parts of the parameter space in the present 
approximation. Note, also, that the effective prediction Y^^^-^'^^f- ( |T8i) is (by definition) a 
positive quantity, therefore the effective description can not reproduce the one-loop result 
at all, which means a large remainder p'^'^--'^<^'^- f |20|) . However, there are a couple of 
circumstances surrounding these situations: first of all, they appear in tiny regions of 
the parameter space; second, and more important, this effect occurs precisely on decay 
channels that have a negligibly small branching ratio, and therefore are phenomenologically 
irrelevant. We can see that the dip in Fig. [TJo around 1250 GeV coincides with the minimum 
of the partial decay width in Fig. [1^. For these reasons we will not try to give a reasonable 
prediction for these decay widths in those corners of the parameter space. From now on 
we will limit ourselves to point out where they appear, so that the reader is warned that 
we can not trust the results in those cases. Outside of this dip, there are two different 
regions. For squark masses larger than 1250 GeV the one-loop correction is around —45% 
whereas the remainder correction (120|) is around —28% that means that, roughly, one third 
of the one-loop corrections can be described as coming from the effective couplings (flGll . 
Since the corrections (from both: one-loop and effective couplings) are quite large, one 
can provide the improved ([21]) description. In the present situation is larger 

than the corrections from the effective couplings and the fixed-order one-loop corrections, 
^Yuk.-imp. ^ —57% but it accounts for two kind of effects: the contribution of higher order 
terms, and the dependence of the radiative corrections on the external momenta. On the 
other side of the plot, for squark masses below 1250 GeV, the situation is quite different. 
The one-loop corrections are relatively small (15%), whereas the effective description gives 
a slightly smaller result, this region has light particles running in the loops, and the one- 
loop functions are expected to depend much on the external momenta. In this situation 
the effective description can not describe properly the radiative corrections. The improved 
description ^M . on the other hand, basically coincides with the fixed order one- loop result. 
In summary: our improved description includes the higher order terms of the Yukawa- 
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Figure 2: a) Top-squark partial decay width into the hghtest chargino and b) relative 
corrections, as a function of the gluino mass. Shown are the different approximations: 
tree-level; one-loop; Yukawa-effective coupling; Yukawa-improved. The input parameters 
are given in eq. ([Ti 



effective couplings, and the external momenta dependence of the one-loop corrections. 
It is able to describe both situations: when the one-loop computation gives a sufficient 
approximation, and when higher order corrections are important and should be taken into 
account. 

Be as it may, it turns out that the effective description using just the Yukawa threshold 
corrections f|T6l) is not enough for the squark decay widths description. The one-loop 
corrections develop a term which grows as the gluino mass rrig [32], which is absent in 
the effective Yukawa couplings ( |T6i) . Fig. [2] shows the comparison of the tree- level, one- 
loop, and Yukawa-effective computations of the top-squark decay width into charginos, as 
a function of the gluino mass, for the input parameters of eq. f|T4|) . One can clearly see 
the log-like behaviour of the one-loop corrections, which can not be reproduced using the 
effective coupling. The logarithmic terms in the gluino mass are a clear example of the 
non-decoupling effects. To understand their origin, one can think on the following: in a 
process where all external (initial and final state) particles belong to the SM sector, one 
can separate the loop contributions in two kinds: SM-like (with only SM and Higgs boson 
particles running inside the loops), and non-SM like (with only SUSY particles running 
inside the loops). Each kind is UV (and IR) finite by itself, so one can remove the non- 
SM-like part from the computation - or take all non-SM-like particle masses very large - 
without any ill effect. In the present case, however, there are SM and SUSY particles as 
external states, the Feynman diagrams can no longer be divided into SM-like and non-SM 
like. The UV divergences from the SM sector are cancelled against UV-divergences of the 
SUSY-sector, and as a consequence if we remove a SUSY particle from the computation, 
the computation is UV-divergent - and therefore meaningless. If we try to remove a 
particle by setting its mass to very large values, this divergence appears as a logarithm of 

■^The variation of the Yukawa-efFective corrections on the gluino mass (fTTl) in Fig. [2^ is of a 0.13%, 
which can not be appreciated in the plot. 
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the corresponding particle mass, and we obtain the aforementioned non-decouphng effects. 
In summary: the QCD corrections to squark decay widths produce exphcit non-decouphng 
terms of the sort logm^. 

To include the logarithmic terms in the effective descriptions, we have extracted from 
the one-loop result the logm^ terms. We have expanded the full one-loop result in the 
limit rrig ^ rriq, and we have used the reduction factors from Ref. [53] to obtain scalar 
quantities. The results for top-squark decays into charginos are: 

A^L'L = f« + rf'L = -fi^« iog!4 . 

ivr /if) 

A^?L = 4*i + *'liL = {A.fiSl'S^ + |4L}log^ . (25) 

here AA^^_ are the full one-loop corrections to the form factors, are the contribu- 
tions from the corresponding one-loop diagrams, SA^_ the counterterm contributions, and 
y4+^_ are the tree-level couplings defined in (fT2|) . fiQ is the scale appearing in the dimen- 
sional reduction of the one-loop UV-divergent integrals, which appears when applying the 
procedure from Ref. [53], and can be though as a renormalization scale. However, these 
expressions do not give any hint at the origin of the logm^ terms, or how can they be 
computed. To understand those terms a renormalization group analysis is in order. 



4 Renormalization Group Analysis 

In order to extract the exact dependence on the renormalization scale, we make a renor- 
malization group analysis, which will allow us to compute the logarithmic terms in rrig. 
To compute those terms, we construct an effective theory below the gluino mass scale, 
which contains only squarks, quarks, charginos, neutralinos and gluons in the light sector 
of the theory, and integrate out the gluino contributions. We find out the renormalization 
group equations (RGE) of the gaugino and higgsino couplings, and perform the matching 
with the full MSSM couplings at the gluino mass scale rrig. In the present computation 
we will consider only logarithmic RGE effects, and neglect the possible threshold effects 
at the gluino mass scale. Since the effective theory does not contain gluinos, only the 
contributions from the gluon have to be taken into account. We will present only the 
computation for the top-squark decay into charginos, the other couplings follow from the 
present computation. 

4.1 Gluon contribution to ^^^j 

The coupling A^^^j is the sum of two terms: a gaugino coupling, and a higgsino coupling, 
as seen in eq. f|T2l) . To shorten up the expressions, we will introduce the shortcuts G, 
which represent the higgsino and gaugino part respectively: 

A^lli = H+ + G ; H+ = -\tRtlyi2 1 G = Ralv*^ . (26) 
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The total gluon contribution to the divergent part of the vertex form factors and the 
wave function renormahzation constants is: 

(Coe/A),^(|.J..ig4L)^^^:L. (27) 

where the first term comes from the vertex form factor's divergent part and the second 
from the fermion and sfermion wave function renormahzation constants. Then, the gluon 
contribution to the (3 function is 

= -2(Coe/A)+ = -^A% = {H^ + G) . (28) 

TT TT 



Therefore, the renormahzation group equation is 

dt TT +'^* 



AlUt) , (29) 



where t = logQ, Q being the renormahzation scale. To solve that equation, we make use 
of the standard RGE for the QCD coupling constant (see e.g. [54]) 

^ - -^^^oiit) . (30) 
Po being the standard QCD /5-function: 

IIN^ - 2Nf Nf 

Po = --^N. , 

3 6 ^ ' 

where iVc = 3 is the number of colors, Nf, Np and Ng are the number of quarks, squarks 
and gluino that have a mass below the scale Q at which we compute the Pq functioE0. 
Inserting expression (|30ll in ( !29l) we obtain 

dAfjt) _ 2 dasjt) 

Solving the equation we obtain 

4Lw _ /«.(^)^* (32) 



Finally, using the QCD running coupling constant 



^ = l-A^,og« . (33) 



"'The RGE evolution is performed by steps, taking into account the change in the /3o function as the 
scale Q crosses thresholds of colored particles. 
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(l-^log^l . (34) 



the running of the right-handed vertex couphng constant is, approximated to 0{as)- 

log TT 

The boundary conditions at Qq = rrig are 

4Lk) = ^+kk)) + gk) , (35) 

then, tlie running couphng constant is 

4L(g) ^ {H^{m,{m,)) + G{m,)) (l - ^ log ^) . (36) 

Now, we have the squark-chargino running couphng constant as a function of the gauge and 
Higgs boson couphngs at the gluino mass scale, but we want to express it as a function of 
couplings at the renormalization scale Q. Note that the gauge part of the coupling G (1261) 
only contains EW gauge couplings, and they do not receive one-loop running contributions 
from the QCD sector, therefore it is a constant term. 

G{mg)=G{Q) = G . (37) 

The higgsino coupling (!26l) . on the other hand, has a dependence on the quark Yukawa 
coupling (or mass), which does run due to QCD corrections, according to the RGE, 

4 

m(g) = m(Qo) f V° , (38) 



\asiQo)/ 

by inserting these expression into eq. fl32|) we can obtain 



4L(g) = A%{m,) 



2 , ,2 



H+[mq[mg)) + G" 



H,(,n,(Q)) I 1 + i^^M log ^ V G (l - log 1 (39) 



vr TTig / \ vr nig 



where, in the last line, we have made the 0{as) approximation. Note that the expressions 
for the higgsino and gaugino couplings are different. Actually, if we write the higgsino and 
gaugino couplings at the scale Q as a function of the couplings at the scale rrig they have 
the same form fl36|) . since they have the same RGE fl32|) . The difference appears when 
we write the higgsino/gaugino couplings at the scale Q as a function of the gauge/Higgs 
couplings at the same scale fl5^ , due to the different running of the gauge fl57|) and Higgs- 
boson fl55]) couplings between the scales rrig and Q. The last line in eq. fl5Ul) agrees with 
the log rrig /fiQ term of the fixed order one-loop expression in (l25ll . 
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4.2 Gluon contribution to A^f^^- 

The gluon contribution to the divergent part of the left-handed couplings A^f^^-, eq. (|T2l) . 
is the same that in the previous case. 



(Co./A)_^^-A™+~^™j^-^™ , (40) 
But now, the coupling only contains a Yukawa like coupling 

= -2(Coe/A)_ = -^ALi = . (41) 

7T 71 



The renormalization group equation is 

dH_{t) _ asit) 
dt TC 

which has as a solution 



H-it) , (42) 

2 

i/_(Q) = i7_(m,(m^))(^#^V° . (43) 



Following the same steps as in the the previous section, we obtain 



log — , (44) 



where in the last line we have made the 0{as) approximation. This expression coincides 
with the higgsino running coupling constant for A^^^- (|39|) . and it also agrees with the 
logrrig/fiQ term of the fixed order one-loop expression in fl25l) . 

4.3 Renormalization Group summary 

The renormalization group running of the coupling constant, can be summarized as follows: 
we can use effective gaugino and higgsino couplings given by, 

9"-{Q) = 9 [ / \ 1 ^log — 



as[mg) J \ 7T rrig 



where \^^^'{Q) are the effective Yukawa couplings defined in f|T6|) . We define then an 
effective partial decay width, computed by replacing g and A in the tree-level expression 
f lT^ by the expressions of eq. 
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Figure 3: a) Partial decay widths and b) relative corrections of the top-squark decay into 
the lightest chargino, as a function of the gluino mass rrig. Shown are different approxi- 
mations to the computation. The effective, improved and remainder computations include 
the logrrig/Q terms. The input parameters are given in eq. ([1 



and we define improved and remainder widths and corrections (r*"*^', S^"^^', S^'^'^'), in 
the same fashion as the Yukawa- effective, Yukawa- improved and Yukawa-remainder of 
eqs. (HHD, dZl- 

The origin of the logm^ terms is a consequence of the SUSY-breaking, and, techni- 
cally, can be seen in different ways, depending on the approximation used to make the 
computation. In the one-loop fixed order computation they appear because of the can- 
cellation of UV-divergences between the gluon and gluino loops (as explained above), but 
in the effective theory point of view, they appear because of the different running of the 
gauge/Higgs boson and gaugino/higgsino couplings. In a fully SUSY theory the gaugino 
(higgsino) couplings are equal to the gauge (Higgs) boson couplings, and they have the 
same RGE, but in a theory with broken SUSY (as the present one) these couplings are no 
longer the same, they have different RGE, and the difference between them is a measure 
of the SUSY-breaking (logm^). 

Fig. [3] shows the partial decay width of the top-squark into the lightest chargino, where 
we include the logm^ terms of eq. (I45ll in the effective description. The input parameters 
are given in eq. dH]). Now the effective description follows the logarithmic behaviour of the 
full one-loop corrections. Moreover, we have checked the validity of our result by comparing 
the full one-loop corrections in the limit rrig S> M with the one-loop expansion of the effec- 
tive description, which contains only log rrig / /ig term^. The results agree with the previous 
results of Ref. [32,34] in the hmit rrig ^ rrig. Our results go beyond the ones of Ref. [32], 
by including the Yukawa terms, including explicitly all the chargino/neutralino and squark 
mixing and couplings, showing the exact dependence from the renormalization group, and 
performing (see below) the numerical comparison with the fixed order computation. 

Fig. [3)d shows the new value of the remainder contributions. At low values of the gluino 

^The comparison is performed by comparing the slope of the different computations in the plots as a 
function of log — . 

^ an 
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mass, the remainder is still large, because rrig ~ rriq and the logarithmic approximation 
does not make sense, but for rrig > 1 TeV the remainder contributions stay below the 2% 
level. 

5 Numerical analysis 

Following the computation and setup of the previous sections, we perform a complete 
numerical analysis. We will concentrate on third generation squark decays (top-squarks 
and bottom-squarks), as they have the richest phenomenology, since only on them the 
higgsino couplings are large enough. The input parameters that we use are given in eq. f|T^ . 
and the resulting spectrum is discussed in section 12. 1[ We have chosen a large value for 
the gluino mass to enhance the effects of the logarithmic terms. 

As a first example we show in Fig. H] the relative corrections to the squark (stop and 
sbottom) partial decay widths into charginos and neutralinos, as a function of the gluino 
mass (rrig). We show the predictions of the effective description including only the effects 
of the effective Yukawa couplings (fT6!) (labelled Yuk.), and including also the log rrig /Q 
terms of eq. (H5!) (lines marked with full circles). We show the partial decays into the 
two charginos and two neutralinos Xi? Xs? the results for the other neutralinos are similar 
to the ones shown. The second neutralino (X2) is mostly a gaugino {w) and its results 
are very similar to Xi^ whereas the fourth neutralino (X4) is mostly of higgsino-type (h) 
and its results are very similar to Xs- plots we observe the same pattern: the 

effective Yukawa couplings f|T6|) do not describe correctly the variation with rrig, only after 
including the log-terms of eq. ( I45l) does the effective description follow the shape of the 
one-loop corrections. The remainder corrections ^"^^"'^•-^em. ( |20|) - those terms that are 
not described by the effective couplings - have a dependence on rrig before including the 
log-terms, but after including them we see that in all channels S^^"^- is essentially flat above 
rrig ~ 1500 GeV, which means that they have absorbed the bulk of the dependence on rrig. 
Moreover, after including the log-terms S^'^"^- is much smaller - between a —2% and a —5% 
- than without them - between a 5% and a —10%. The effects of the log-terms are more 
visible in the gaugino-like channels, where the Yukawa couplings play no role, and the bulk 
of the corrections corresponds to the log-terms. In these channels the corrections change in 
a range of a 10% in the interval rrig = [500, 5000] GeV. In the higgsino-like channels their 
importance is less apparent. On one side the effective Yukawa couplings carry the bulk of 
the corrections (around a —30%), and on the other side the sign difference among the two 
contributions ( l45l) make them to partially compensate each other. But also in this case the 
effective Yukawa couplings alone ( fT6|l do not describe correctly the m^-dependence, and 
there is a range of a 3% variation in the corrections in the studied m^-interval. The spikes 
that are seen in some plots correspond to the threshold singularities for the opening of the 
squark decay into gluinos {q —>■ gq), we recall that in this work we are not interested in 
the region where the gluino decay channel is open {rrig + rriq < rriq), since in that region 
the strong decay is the leading one, and the chargino/neutralino channels have a negligible 
branching ratio, and are phenomenologically irrelevant. Nevertheless we include the plots 
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Table 1: Value of the remainder corrections ( |20l) for specific squark partial decay widths, 
including only the effective Yukawa couplings, eq. (fT6l) . and results including the logrrig/Q 
terms, eq. (H5l) . for tan/5 = 50, and the rest of input parameters as given in eq. (HM . 



also in that region to show the trend of the corrections. 

Next, we show in Fig. [5] the evolution of the different corrections as a function of tan (3. 
Here we see large negative corrections, growing with tan/3. The origin of the negative 
corrections is twofold: on one side the standard QCD running of the quark mass reduces 
significantly the Yukawa coupling, and on the other, for positive values of the higgsino 
mass parameter fi the contributions to Arrig f|T7j) are also positive, decreasing even more 
the effective Yukawa couplings f|T6|) . However, even after taking into account these two 
sources of corrections, still there is a large remainder of relative corrections (up to —67% 
at tan/5 = 50 for some channels) which can not be accounted for. After including the 
logarithmic terms, eq. ( H5l) . the situation is quite different. Now, the effective description 
can reproduce quite well the one-loop results, and the remainder corrections (those that 
can not be described by the effective couplings) are reduced at a level (in absolute value) 
below the 20%. Table [1] shows the value of the remainder corrections (!20l) for all studied 
channels at tan/3 = 50, including the full effective description, and including only the 
Yukawa corrections. We see that in all channels the corrections are reduced significantly 
after including the logarithmic terms. 

Finally Fig. [6] shows the evolution of the corrections to the partial decay widths as a 
function of the SU{2) l squark mass parameter Mj^. The abrupt change which is seen at the 
middle of the plots corresponds with the situation in which Mj^ ~ ^/r' physical 
states suffer an abrupt change between left and right chirality. This explains the difference 
in value and behaviour of the corrections in the regions of Mj^ below and above that point. 
This is also the situation, shown in Fig. [1], where the partial decay widths can become 
zero, and the one-loop corrections can become non-perturbative. In all situations the 
description including the log-terms provides a better description of the radiative corrections 
with a S^^"^' f l20p much smaller than with the effective Yukawa couplings (1161) alone, and 
a much softer variation, meaning that the description of eq. f H5|) is accurate for all values 
of the squark mass. Let us remember, that by changing the SUSY parameter Mj^, the 
physical squark masses also change, and that since the renormalization scale is taken to 
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be the decaying squark mass, the log-terms of eq. (H5|l effectively run also with Mj^ as 
~ logrrig/Mj;^ (for mostly- left-handed squarks) - a contribution that can not be described 
with the effective Yukawa couplings (fT6ll . 

6 Summary and conclusions 

We have proposed and analyzed an effective description of squark interactions with charginos 
and neutralinos in the MSSM. We have applied it to the partial decay widths of squarks 
into charginos and neutralinos. We have compared it with the full one-loop corrections, 
and have proposed a way to combine the effective description (which includes higher order 
terms) with the complete one-loop description (which includes all kinetic and mass-effects 
factors), providing an improved computation, r*'"^' f|T9|) . The difference between the effec- 
tive description and the improved computation is encoded in the remainder contribution, 
^rem. f[2Ul) . which gives a measure of the precision of the effective description. 

The effective description includes the effective Yukawa couplings flTBl) . which take into 
account the resummation effects [49-52]. Note that the computation of the threshold 
corrections (fTTI) includes only the (tan/5, cot/3) proportional terms, since the terms that 
would be proportional to the trilinear couplings {Ab, At) are actually subleading [52]. This 
description produces large remainder §Yuk.-rem. corrections, and does not reproduce the 
behaviour of the one-loop corrections as a function of several parameters - notably it is 
missing a logm^ term. Therefore, we conclude that it does not reproduce satisfactorily the 
one-loop corrections, and it is not a good approximation. 

We have computed the missing log rrig terms with the help of the renormalization group 
- eq. (jlS]) and found agreement between the renormalization group analysis and the 
large mass expansion of the one- loop result fl2S]) . After including also the logm^-terms of 
eq. ( H5|) . the effective description produces a reasonable approximation to the radiative- 
corrected partial decay widths of squarks into charginos and neutralinos, as shown by a 
small absolute value of the remainder contributions S^'^"^', and by a nearly-flat behaviour 
of the corrections as a function of different parameters. 

The origin of the logarithmic terms can be explained in different (complementary) ways, 
depending on the kind of approximation that we take. First of all, from the fundamental 
point of view, they are non- decoupling terms that appear due to the supersymmetry- 
breaking. Since we are testing SUSY relations (equality of the gauge/ Yukawa couplings 
to the gaugino/higgsino couplings), and SUSY is broken (by the term rrig among others), 
we have to find some effect that tells us about the breaking of SUSY at that scale - e.g. 
a logm^-term. Second, from a fixed-order one-loop description point of view, the SUSY 
relation appears because the UV-divergences of the loops containing gluinos cancel with 
the UV-divergences of the loops containing gluons, and the log-terms that accompany 
those loops combine between them - producing a logm^-term. And third, from a renor- 
malization group (and effective theory) point of view, the gaugino/higgsino couplings run 
different than the gauge/Yukawa couplings in the region where SUSY is broken, that is, 
for scales below rrig. The running of the gaugino/higgsino couplings from the scale rrig to 
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the chosen renormahzation scale produces logmg/Q-terms. Since the Yukawa couphngs 
already include some QCD running, whereas the EW-gauge couplings do not, the relation 
between the gauge/gaugino and Yukawa/ higgsino couplings has some differences ( l45l) . 

The presence of the non-decoupling logm^-terms implies a deviation of the equality 
between the higgsino/gaugino and Higgs/gauge couplings predicted by exact SUSY. This 
deviation is important, and has to be taken into account in the experimental measurement 
of SUSY relations. At the same time, it gives us access to information about heavy particles 
that can not be directly produced at the LHC/ILC. For these reasons it is important to 
include these effects in the computation for the predictions of squark observables at the 
LHC and the ILC. The effective description of squark/chargino/neutralino couplings given 
by eqs. f|T6|) . (fT7|) . (145!) . is simple to write, and to introduce in computer codes, it costs 
little computational power, and provides a reasonable description for squark decays into 
charginos and neutralinos, so it can be used in monte-carlo generators and other computer 
programs that provide predictions for the LHC and the ILC to improve their accuracy at 
a minimum cost. 
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